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1. Introduction 
We study the operator L pcT defined by 
(1.1) L pcl u := ie(f (x)u (x)) 1 + iu'(x) 

in which / is a given 27r-periodic function having the following properties: 
(1-2) f(x + n) = -f(x), f(-x) = -f(x); 

and also 

(1.3) f(x)>0 forxe(0,7r). 

In particular it follows that /(7rZ) = 0. We assume that / is continuous, and 
differentiable except possibly at a finite number of points, the points of non- 
differentiability excluding irZ. We assume that /'(0) = 2/tt and that < e < it. 
We consider (jl.ljl on the domain 

(1.4) V = {u G i 2 (-7r,7r) | LperW G L 2 (-ir, w); u(-ir) = u(ir)}. 

Remark 1. Of course it is not obvious that functions u G L 2 (—ir,ir) such that 
L pm u £ L 2 {— 7r,7r) have boundary values u(±tt); this was proved in [BLMj . where 
we showed that if u G L 2 (—ir,ir) and L pcr u G L 2 (~ir,ir) then u G H 1 ^— ir, w). 

The main results of this paper are 

Theorem 2. TTie spectrum of operator L pcl is 

(a) rea/ 

(b) purely discrete, i. e. it consists only of isolated eigenvalues of finite multi- 
plicity with no accumulation points apart from, possibly, infinity. 

Part (a) has been partially proved in [BLMj . where we showed that all the 
eigenvalues (if they exist) are real. The rest of Theorem [2] follows from 

Theorem 3. The resolvent (L pcl — A) -1 is a compact operator on L 2 (— tt,tt) if A 
is not an eigenvalue of L peT . 

Remark 4. The spectrum is always non-empty, as zero is an eigenvalue correspond- 
ing to a constant eigenfunction. 
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In order to prove Theorem [3j we show that when A is not an eigenvalue of L pe r 
(which is guaranteed, for instance, if A is not real) then the boundary value problem 

(1.5) ie(f(x)u (x)) + iu (x) — \u(x) = F(x) — ir < x < ir 

with periodic boundary conditions u(— ir) = u(ir) has a unique solution u G T> for 
every F £ L 2 (~ir, it). 

The compactness of the resolvent is demonstrated by an "explicit" construction 
of a bounded Green function G(x, s) such that u(x) = J_ G(x, s)F(s) ds. The 
properties of G are established by studying the solutions of an associated homoge- 
neous equation in Sections [2] and [3] 

Motivation and scope of the present paper. Our interest in the operator (|l.ip 
with domain determined by condition (|1.2jl is primarily motivated by BcO'BSa 
and [Da2 , where f{x) — (2/7r)sinic, and therefore (|1.5p takes the form 

(1.6) ie (sm(x)u'(x))' + iu'{x) — \u{x) — F(x), 

with < e < 2. This equation arises in fluid dynamics, and describes small oscil- 
lations of a thin layer of fluid inside a rotating cylinder. From a purely theoretical 
perspective, the eigenvalue problem associated to L per has recently drawn a sub- 
stantial amount of attention: see |ChPe] . |Welj . |We2] . |DaWej and |Ch KaPy| . 
Despite the fact that (|1.6I) is highly non-self-adjoint, the spectrum of (|1.6p consists 
exclusively of real eigenvalues of finite multiplicity, it is symmetric with respect to 
the origin and it accumulates at ±oo. It is also known that the eigenfunctions do 
not form an unconditional basis of L 2 (—tt,tt). Moreover, Davies and Weir [DaWej 
have recently found explicit asymptotics of the eigenvalues as e — » 0. 

In |BLMj we established that the eigenvalues of L per are all real and form a 
symmetric set with respect to the origin. Theorem [5] above rules out completely 
the possibility of a non-empty essential spectrum and it answers an open question 
posed in our previous paper. 



2. A NAIVE FROBENIUS ANALYSIS 

Let p satisfy the integrating factor equation 
Then u(x) is a solution of (|1.5p iff 

(2.2) {p{x)u'(x))' + l ^ju{x) = - (*) - t < x < TT. 

For future use, we also recall the homogeneous differential equation 

(O) (p(x)u'(x)Y + —^-U(x) = -TT<X<TT 

In order to understand how p behaves near x — and x = n it is useful to consider a 
simple model. Suppose that near x = 0, the function / satisfies f(x) — 2x/tt - this 
maintains the normalization /'(0) = 2/ir. Then (|2.ip yields logp = log(2x/7r) + 
7r/(2e) log(x), whence p(x) — Cx 1+C ^ where C is an arbitrary non-zero constant 
and c = ir/2. Similarly, near x = ir, we can consider a simple model f(x) = 
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(2/tt)(tt — x) and obtain p(x) — C(tt — x) 1 ~ c l' L . In |BLM] we proved that, under 

the minimal assumptions (|1.2p on /, the following results hold: 

(2.3) 

, . / x 1+c ^ x ~ p(x) / x c l c x ~ 

PW ~ \ (tt - x) 1 "^ x /(xj ~ \ (tt - x)- c A x ^ tt. 

By considering the model f(x) = 2x/ir, p(x) = x 1+c / c in a neighbourhood of 
the origin and looking for solutions of the differential equation (|2.2I[) in the form 
u(x) — x u {\ + a\x + a-ix 1 + • • • ) one can establish the asymptotic behaviour of 
solutions for this model in a neighbourhood of the origin; similarly near x — tt. In 
BLM we show that, under hypotheses (|1.2I) on /, there exist solutions u of (|2.2ip 
such that 

f x" c / e or 1 x ~ 

«(x, A) ~ < 1 or (tt — x) c / c x ~ 7r 
[ 1 or (i + 7r) c / e X ~ — 7T 

This implies the existence of a unique solution </>(x) = (f>(x, A) of (12.2ip . such that 

(2.4) (j>(x,\)~l, x^O. 

The following result is crucial for reducing the problem to the interval (0, 7r); in a 
sense it plays the same role in the analysis as the orthogonal splitting of the infinite 
matrix operator in Davies [Da2j in £ 2 (Z) into three operators, hence reducing the 
problem to a problem in £ 2 (N). 

Lemma 5. The solution 4> has the symmetry property 

(2.5) 0(-x,A)=#s,-A). 

Proof. Define a function v(x) = 0(— x, A). A direct calculation shows, thanks to 
the symmetry conditions (|1.2[) . that v satisfies (|2.2lj) but with A on the right hand 
side replaced by —A. It also satisfies v(0) = 1. However eqn. (|2.2ip has only one 
solution with this property, namely 4>{x, —A). This proves the result. □ 

We emphasize the important fact that A is an eigenvalue of (|2.2ip if and only if 
<p possesses the additional symmetry property 

(2.6) <K-7r,A)=0(7r,A). 

In [BLM] we also show that there is also a second solution tp(x) = %[)(x, A) of 
([2~2f]) satisfying 



(2.7) ^(x,A) 



|x|~ c / £ x~0 

|X =F 7r| C / £ X ~ ±7T 



Observe that ^>(— 7r, A) = ^(tt, A) = 0, and that VK^i A) blows up when x ~ 0, at 
least when A is not an eigenvalue. 

Consequently, when A is not an eigenvalue, we can also normalize ipix, A) by the 
condition 

p(x)tp' (x)<f)(x) — p(x)(j>' (x)ip(x) = 1 — tt < x < TT. 

The Wronskian in the right-hand side here is obviously a constant, and below we 
will always assume that ip(x, A) satisfies this condition. 
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3. Construction of the Green function 



Now, we are back to constructing explicitly the L 2 solution u(x, A) of (|2.2j) 
assuming that (|2.6[) fails. By the standard variation of parameters technique the 
general solution of (|2.2I) takes the form 



u(x, A) = ip(x, A) 



-lp 



s,X)[-^F)(s)ds 



(3-1) 



-<■>{.,: \) I I '(.s. A) ( -^F ) - 



+ A<£(x) + B^(x) , 

where A and B are arbitrary complex constants. 

It remains to check that one can choose constants A and B in such a way that 

(a) u e L 2 (-ir, 7r), 

(b) ip 0r u £ L 2 (—TT 7 ir), and 

(c) it(— 7T, A) = w(7r, A). 

It is in fact sufficient, and easier, to check that one can choose constants A and B 
such that 

(a') u(x, A) is continuous at x = 0, 
(b') ie/u' + iu is continuous at x = 0, 

and (c) all hold. Indeed, (|3.1j) together with (a') implies (a), and together with (b') 
and (c) implies (b). 

Note first, that by (a') and the behaviour of ip near the origin, one is tempted 
to take B = in (|3.1[) . We shall show that this choice is indeed the right one by 
the careful analysis of the remaining terms in (|3.1|1 . 

By the Cauchy-Schwarz inequality and (|2.3j) . (|2.4|) . (|2.7p we have, 



c(.r.A) / ,,(.s.A)(-^F] (.s)d.s 



(3.2) 



< \ip(xA)\ 



\4>{s,x)\ 



-ip 



< C\x\ 

< C|: 



-c/e|„|c/e+l/2 



[ X \F(s)\ 2 d S 
Jo 



2 \ V2 

ds 



|F(s)| 2 ds 



1/2 



1/2 



|!/2| 



Here, and throughout the rest of this paper, C denotes a generic positive constant; 
|| • || is the standard norm in L 2 (—ir, tt). 
Similarly, 



<»•<•• A. / r(.s.A) ( -^-F) (s)d.s 



(3.3) 



< |^,A)| 



|V(s,A)| 



1/2 



ds 



\F(s)\ 2 ds 



1/2 



< C- 1- \n-x\ 1/2 
<C\n-xnF\\ 



£\F( S )\ 2 ds 



1/2 
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so it is bounded. It is also continuous as the integrand is a product of a bounded 
function ipp/ f and an L 2 function F. 

Also, 4> is continuous at zero, and ip is not, and therefore (a') holds if and only 
if B = 0. 

To check the condition (b'), it is now sufficient to verify that fu' is continuous 
at zero. Differentiating (|3.1|) with respect to x gives 



/"'(.•''■ A) : ■ /<•'(,-. A) / o(.s.A) ( -^F ] • 

+ A) / V(s,A) ( ^F) (s)ds 



(3.4) 



W(a:) , 



as the contributions from differentiating the integrals cancel out. 

The last two terms go to zero as x — > 0, and we only need to check the continuity 
of the first term. Similarly to Q3.2p . we get 



Hr)c'(r.X) j „(.s.A) ( y^/' S 



(3.5) 



< \f(xW{x,\)\ 



|0( S ,A)| 2 



-ip 



ds 



< C\x\\x\- c/e - l \x\ c/e+1/2 



<C\x\V*\\F\\, 



\F{s)\ 2 ds 



1/2 



1/2 



|F( s )| 2 d S 



1/2 



which proves (b'). 

Finally, we need to guarantee that we can choose a value of constant A to ensure 
that condition (c) holds. Direct substitution, again taking account of (12. 3[) . (|2.4|) . 
([2~7]) gives 



u(tt, A) = A4>(tt, A) , 



-ip 



u(-7r, A) = A^-tt, A) + 0(-tt, A) y ^(s, A) [-jfF) (*) ^, 

and as A is assumed not to be an eigenvalue, and therefore (|2.6[) is not satisfied, we 
can choose 



(3.6) 



.4 



<K*,A) 



-iff 



F (s) ds 



i(7T, A) 
H^A) 



f 



This proves the existence of the resolvent of (L per — A) -1 : L2{— ir, ir) — > 7r, n) 
for each A which is not an eigenvalue of L pel . Thus, the spectrum of L per is pure 
point and real. 

Now we proceed to writing down the expression for the Green function G(x, s) 
of Lper- We note that (|3.ip can be written, with account of B — and (|3.6p . as 



u(x,X)= G(x,s)F(s)ds = (Gi(x,s) + Gn(x,s) + Giu(x,s))F(s)ds 
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where we set G(x, s) := Gi(x, s) + Gn(x, s) + Gm(x, s) , with 
>(^A)^( s ,A)(^gi) ifM>N, 



Gi(x, s) 
G u (x,s) 



otherwise; 
otherwise; 



Gni(x, s) := ^(as, A)^(s, A) ( -^f) / ( 4r^C\ ~ 1 ) for all - ^ < s < tt . 



The functions Gu(x,s) and Gm(x,s) are bounded by (|2.3p , (|2.4[) . (|2.7p , so we 
need to look at Gi(x, s). The only scope for trouble in the expression for G\ lies in 
the fact that ip(x, A) blows up as x — > 0. However if x is small then, in the region 
\x\ > \s\ where G\ is nonzero, (|2.7p and (|2.3p yield 

< C\x\- c/e \s\ c / e < C. 

Thus G is bounded and hence is the kernel of a compact operator on L 2 (—ir, n). 
Thus Lper has compact resolvent and purely discrete real spectrum. 



p(s) 



4. SCHATTEN CLASS PROPERTIES OF THE GREEN FUNCTION 

Wc first recall the standard notion of Schattcn class operator. Let T be a compact 
operator and consider its (oo-dimensional) singular value decomposition (see e.g. 
[GoKr]): 

OO 

T = Y, a j\vj)(wj\ 

3=0 

where the singular values ctj > and the two sets of vector {vj} and {wj} are 
orthonormal and not necessarily equal. Here we use the bra-ket notation: \v) (w\u = 
(u, w)v. For p > 0, we say that T is in the p-Schatten class, T £ C p , if 

\\T\\ p :=(^£a p )j 

Note that C\ are the trace class operators and C2 are the Hilbert-Schmidt operators. 

Theorem 6. The resolvent (A — L per )~ l is in C p for all p > 1. 

Proof. Let G\(x, s), G\\(x, s), G\u(x, s) be the components of the Green function 
G(x,s). It suffices to show that each of the corresponding integral operators 
i?i(A),i?„(A),i?in(A) isinC p . 

We start by showing that Rm(X) € C±. Indeed, Gm(x, s) — v(x)w(s) where 

A) and w {s)=^s,X)^^j/(^^L- 

By virtue of (O, and both v,w G L 2 (-7r,7r). Then i? m (A) = \v)(w\ 

and ||i?in(A)||i = ||-u|| i 2 { _ 7r:7r) |MU 3 (-7r,7r) < 00, as required. 
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Let us now show that -Rn(A) £ C p for all p > 1. Let fin = {(x, s) € [— tt,tt] 2 : 
x < s}, be the supports of Gu(x, s). Decompose the characteristic function of flu 
as 

oo 2 j -l 
0=0 i=0 

where there are 2 J+1 intervals 
i i+l 



2i- 



it: (i + 1)tt 

' 7T) — 7T 



2J+ 1 



, z = 0,...,2^ +1 -l. 



Then 



oo 2 j -l 

G n (x,s) = t nu (x, s)v(x)w(s) = E E v(x)t l2ij (x)w(s)l l2i+lj (s). 

j=0 i=0 



Let Si 



\vt l2i J(wti 



Then 



■= \\Si,j\\ = IMIl2(/ 2i , 3 )IMIl2(/ 2i+1iJ ) 



< 



23 



where m is independent of i and j. The constant m depends on A and it is finite 
as a consequence of (|2.4[) . ()2.7[) and ()2.3|) . Let 



(4.1) 



Si 



E = E 



i=0 



i=0 



vl 



wt 



l«lli 2 (/2i,j)/ \ IFIIi 2 (Wi,j) 

Since the intervals ijj are pairwise disjoint for a fixed j, the right side of (|4.1|) is a 
singular value decomposition for Sj. Then 

i/p 

m 



<2 D -1 



J 3\\P 



< 



i=0 



p-l \ 3 

2~ 



By the triangle inequality, this ensures that 

||i?ii(A)|| p < oo 

for all p > 1 as required. 

The fact that i?i(A) £ C p for all p > 1 follows by an analogous decomposition of 
the support Q\ = {(x,s) £ [— tt,tt} 2 : \x\ > \s\} as the union of disjoint rectangles 
and a very similar argument. 

□ 

Remark 7. Let A„ be the eigenvalues of L per and A ^ A„. By virtue of [DuScl 
cor.XI.9.7], the series J2^Lq(^ — A„)~ p converges absolutely and 

OO 

(4-2) ElA-A^-^IKA-Lper)- 1 !!? 

3=0 

for all p > 1. According to the results of |We2) on the case f(x) — (2/tt) sinx, it 
is known that A„ ~ n 2 as n — > oo. Hence we know that (A — L per ) _1 Ci/2- As 
L pcr ^ L* el , the inequality in (j4.2|) can not generally be reverse for any p > 0. The 
question of whether (A — L pcr ) _1 £ C p for p > 1/2 will be addressed in subsequent 
work. 
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